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Abstract 

We study the global flow of the anisotropic Manev problem, which 
describes the planar motion of two bodies under the influence of an 
anisotropic Newtonian potential with a relativistic correction term. 
We flrst flnd all the heteroclinic orbits between equilibrium solutions. 
Then we generalize the Poincare-Melnikov method and use it to prove 
the existence of infinitely many transversal homoclinic orbits. Invok- 
ing a variational principle and the symmetries of the system, we finally 
detect infinitely many classes of periodic solutions. 



1 Introduction 

The anisotropic Manev problem describes the motion of two point masses in 
an anisotropic configuration plane under the influence of a Newtonian force- 
law with a relativistic correction term. The isotropic case is the classical 
Manev problem; its origins lie in the work of Newton, who introduced it in 
Principia aiming to understand the apsidal motion of the moon (see [|ll|, 0]). 
Manev found in the 1930s that a proper choice of the constants that show 
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up in the correction term allow the theoretical explanation of the perihelion 
advance of Mercury and of the other inner planets. 

The first author suggested the study of the anisotropic Manev problem in 
1995, hoping to find connections between classical, quantum, and relativistic 
mechanics. It was indeed proved in [jlO| that the rich collision-orbit manifold 
of the system exhibits classical, quantum, and relativistic properties. This 
encouraged further studies, as for example [0 and In [jl5|, using a suit- 
able generalization of the Poincare-Melnikov method (see 0, ^ ^ for 
the classical approach or ||, |[ for a parallel, at least in part, complementary 
approach), we proved that chaos occurs on the zero-energy manifold, thus 
showing the complexity of the dynamics. Using perturbations techniques and 
the Poincare continuation method, the second author investigated in the 
classes of periodic solutions that arise from symmetries in the case of small 
values of the anisotropy parameter. 

In this paper we gain a better understanding of the complicated global dy- 
namics encountered in this problem. We first prove that negative-energy 
solutions are bounded and find the heteroclinic orbits that connect the equi- 
libria of the collision manifold, which we obtain through McGehee-type trans- 
formations (see |2^). Physically they correspond to ejection-collision orbits. 
Then we employ perturbation techniques to detect possible global chaotic 
behaviour. As remarked in |2^, the perturbation analysis of []T5|, ^ cannot 



be used to study ejection-collision solutions. However, we surpass this diffi- 
culty with the help of McGehee-type coordinates, which allow us to view the 
anisotropic Manev problem as a perturbation of the classical Manev case. 
Using an approach inspired by P] , which works in some degenerate cases — 
as for example those of unstable nonhyperbolic points or critical points lo- 
cated at infinity (see [0, ||, ||, |l^), we develop a suitable extension of the 
Poincare-Melnikov method, which we use to prove the existence of transver- 
sal homoclinic orbits to a periodic one. It is interesting to note that our result 
extends the one obtained in |^, ^ to a non-Hamiltonian system that has 
negatively and positively asymptotic sets to a nonhyperbolic periodic orbit. 
In the present context the asymptotic sets are the stable and the unstable 
manifolds. 

Then we return to the original coordinates and apply a variational principle 
for detecting periodic orbits. Using the rotation index, we divide the set of 
periodic paths into homotopy classes, which are Sobolev spaces. Then we use 
the lower semicontinuity version of Hilbert's direct method (due to Tonelli, 
see [26 1) to find a minimizer of the action in each class. According to the 
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least action principle, the minimizer is a solution of the anisotropic Manev 
problem. We prove that the minimizer exists, belongs to the homotopy class, 
and is a solution in the classical sense. This generalizes a result obtained by 



the second author, [^, where it was shown that such orbits exist for small 
values of /i > 1. In the end we put into the evidence some new properties of 
symmetric periodic orbits. 

The idea of using variational principles to obtain periodic orbits for n-body- 
type particle systems first appeared in and has been recently used in 
connection with symmetry conditions to obtain new periodic orbits in the 
classical n-body problem (see 0). But unlike the Newtonian case, the Manev 
force is "strong" (as defined in so the variational method is easier to 

apply in our situation than in the Newtonian one. This is because in the 
Manev case we do not have to deal with the difficulty of avoiding collision 
orbits, which have infinite action and therefore cannot be minimizers. 
Our paper is organized as follows. In Section 2 we write the equations of mo- 
tion and transform them to an equivalent system using a "blow-up" technique 
devised by McGehee, which allow us to introduce the concept of a collision 
manifold. In Section 3 we present two global results: the boundedness of the 
solutions for negative energy and the existence of certain symmetric ejection- 
collision orbits. In Section 4 we describe the anisotropic Manev problem as a 
perturbation of the Manev case. In Section 5 we develop a suitable general- 
ization of the Poincare-Melnikov method and in Section 6 we apply it to find 
infinitely many transverse homoclinic orbits that show that the dynamics of 
the problem is extremely complex, possibly chaotic. Finally, in Section 7 we 
use a variational principle to prove the existence of infinitely many classes of 
symmetric periodic orbits. 



2 Equations of Motion 

The (planar) anisotropic Manev problem is described by the Hamiltonian 

rr^l b__ 

where /i > 1 is a constant, q = (x, y) is the position of one body with 
respect to the other considered fixed at the origin of the coordinate system, 
and p = {px,Py) is the momentum of the moving particle. The constant /i 
measures the strength of the anisotropy and we can very well take yU < 1; 
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but to remain consistent with the choice made in previous papers, we will 
consider /i > 1. For /i = 1 we recover the classical Manev problem. The 
equations of motion are 

The Hamiltonian provides the first integral 

H{p{t),qit)) = h, (3) 

where /i is a real constant. Unlike in the classical Manev case, the angular 
momentum K{t) = p(t) x q(t) does not yield a first integral. This is because 
the anisotropy of the plane destroys the rotational invariance. 
Since our first goal is to study collision and near collision solutions, it is 
helpful to transform system (0) using a method developed by McGehee pp 



The idea is to "blow-up" the collision singularity, replace it with a so-called 
collision manifold and extend the phase space to it. The collision manifold is 
fictitious in the sense that it has no physical meaning. However, studying the 
flow on it provides useful information about near-collision orbits. Consider 
the coordinate transformations 

r = |q| 

9 = arctan(|//a;) , . 

V = rr = (xpx + yPy) 
u = r'^9 = {xpy - ypo,), 

and the rescaling of time 

dr = r'^dt. (5) 

Composing these transformations, which are analytic diffeomorphisms in 
their respective domains, system becomes 

r' = rv 
v' = 2r^h + 
9' = u 

u' = (l/2)(/i- l)(rA-3/2 + 26A-2)sin2^ 
and the energy relation takes the form 

u^ + v^- 2r A-^/2 _ = 2r^h, (7) 
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where A = /i cos^ 9 + sin^ 9 and the new variables (r, f , 9^ u) G (0, oo) x M x 
S"^ X M depend on the fictitious time r. The prime denotes differentiation 
with respect to r. 
The set 

C = {(r, V, 9, u)\r = and the energy relation (0) holds} (8) 

is the collision manifold, which replaces the set of singularities {(q, p) |q = 0}. 
This 2-dimensional manifold, embedded in M.^ x S^, is homeomorphic to a 
torus and it is given by the equations 

r = and + = 2bA~\ (9) 



The flow on the collision manifold was studied in detail in [|TU]. Here we 
will briefly recall its main features. Let's consider the restriction of sys- 
tem @ to C. The solutions of the restriction lie on the level curves v = 
constant of the torus C. There are eight equilibrium points. In the vari- 
ables {r,v,9,u) the first four equilibria are = (0, ±a/26/ /i, 0, 0) and 
= (0, ±A/26//i, TT, 0). The corresponding eigenvalues are real and take the 
values ±A/26//i, 0, ±a/26(1 — The other four equilibria are ^^^/2 ~ 



(0, ±\/2fe, ±7r/2, 0) and the corresponding eigenvalues are 
where the last two eigenvalues are purely imaginary since /i > 1. Moreover 
there are eight heteroclinic orbits which lie in the level sets v = ±^26/ /i. 
All the other solutions are periodic (see Fig. 1). 



3 Heteroclinic Orbits and Bounded Solutions 



The flow near the collision manifold was studied in |T2|, in which most of the 
results are essentially local. In this section we will prove two global results 
that extend the understanding of the problem under discussion. The first 
one concerns the boundedness of solutions on negative energy levels. 

Theorem 1 For any negative value of the energy constant, h < 0, there ex- 
ists a positive real number M such that any given solution (r(r), t;(r), 9{t),u{t)) 
of system satisfies the relation r < M. 

Proof: Let us assume that there is no M with the above property. Then 
at least one unbounded solution exists. Since by the energy relation 
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Figure 1: The flow on the colhsion manifold, which is formed by periodic 
orbits, eight equilibria, and eight heteroclinic orbits. 

+ i;^ = 2r^/i + 2rA^^/^ + 26A^^, and since h < 0, there is some r = r(r) 
such that + f ^ is negative — a contradiction. This completes the proof. 

The next result deals with the existence of heteroclinic orbits connecting the 
equilibria but lying outside the collision manifold. But before stating and 
proving it, let us recall some facts that summarize the behavior of the flow 
near the collision manifold. Denote by Prj the periodic orbit on C having 
V = T]. The following property was proved in |]TU[. 

Proposition 1 On the collision manifold C the equilibria and are 
saddles whereas the equilibria A^^j^ are centers. Outside the collision mani- 
fold the equilibria A^, A~^^^^, and A'^ have a 1-dimensional unstable analytic 
manifold, whereas the equilibria A^ , A±t^i2, and A~ have a \- dimensional sta- 
ble analytic manifold. Each periodic orbit on C with v = ri>0(v = 'r]<0) 
has a 2-dimensional local unstable analytic manifold, while the periodic orbit 
= has both a 2-dimensional local unstable and a 2-dimensional local stable 
manifold (see Fig. 2) 

The above properties are local, the following one, however, is global. We will 
now show that the equilibria with positive v coordinate have a 1-dimensional 
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global unstable manifold while the equilibria with a negative v have a 1- 
dimensional stable manifold. Moreover, the equilibria are connected by het- 
eroclinic orbits starting from an equilibrium with positive v and ending in 
the symmetric one with respect to the {6, u) plane. 

Theorem 2 There are four heteroclinic orbits outside the collision manifold 
C : 7-7r/2 ;7o n-K/2 ,1tt Connecting respectively A^^i^ with ^1^/2' '"^^^^ ' 
A^l^ with A~^^, and with A~ (see Fig. 2) 

Proof: First we show that u = and 6 = 0, vr, ±7r/2 describe four invariant 
sets. Consider 6{0) = 6o = and n(0) = uq = 0, as initial conditions. Then 
9 = 0, u = satisfies system (^), hence (by the uniqueness property for 
solutions) u = 0, 6 = define an invariant set. The same reasoning can be 
applied if ^ = ±7r/2 or ir. 

Now let's study the energy relation (|^) when n = and 6* = 0, tt. After 
simple computations we get 

+ ( ^/2\h\r - -7^] =- + An- (10) 

The above equation describes an ellipse whose intersections with r = give 
V = ±-^26/ /i, which are exactly the equilibrium points and A^. Moreover 
the maximum value of If I is 

(11) 

and the maximum value of r, attained when f = 0, is 











2\h\J 



1(1 -Ahb) 



''""'^ ~ i2^\h\) + 2\h\ ■ ^^^^ 

Consequently for m = 0, 6' = (6^ = tt) there exist heteroclinic orbits 70, (7,r) 
ejecting from A^ {A^) and tending to A^ (A^) (see Fig. 2). 
Similarly when n = and 9 = ±7r/2 the energy relation can be reduced to 
the form 
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which describes an elUpse. The intersections with r = are v = and 
represent the equihbria A^^^^. In this case 

Vma. = V26+l/2|/l| (14) 

and 



1 , ^^^^ 



2\h\ 2\h\ ■ 

Thus we found heterochnic orbits 7±7r/2 ejecting from A^^^^ ^^'^ tending to 
^±7r/2 (^^^ -^^S- 2)- '^^^^ completes the proof. 




Figure 2: The flow can reach the collision manifold at the equilibria or at 
any of the periodic orbits. There are four heterochnic orbits 7_7r/2 ,70 577r/2 
,7^ connecting respectively A~^^^^ with with Aq , A^^^ with A'^^, 

and At with Az . 



4 A Perturbative Approach 

We will now write the anisotropic Manev problem as a perturbation of the 
classical Manev case. Consider weak anisotropics, i.e., choose the parameter 
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/i close to 1. Introducing the notation /i — 1 = e > with e ^ 1, we can 
expand the equation of motion in powers of e to obtain 



r' = rv 

v' = 2r^/i + r — e(r/2 cos^ 



9' = . 

u' = e(r/2 + 26) sin 2^. 

The energy relation becomes 

u"^ + -2r -2b + e(r + 2b) cos^ 6 = 2r^h. (17) 

For e = 0, system (|16D and equation (0) yield the Manev problem. The 
collision manifold is the set of solutions given by 

r = 0, u'^ + v^ = 2b. (18) 

Notice that, from the geometric point of view, the collision manifold is a 
cylinder in the three-dimensional space of coordinates {u, 6, v) and, since 
6 G [0,2-7r], it follows that this cylinder can be identified with a torus. The 
flow on the collision manifold is formed almost exclusively by non-hyperbolic 
periodic orbits, except for the upper and lower circles of the torus given by 
r = 0, u = 0, V = ±^/2b, which consist of equilibrium points. There is 
only a single orbit ejecting from each fixed point of the upper circle v = \/26 
and a single orbit tending to the lower circle v = —\/2b (see [|TB|). Moreover 
it can be easily proved (see |T^) that for every periodic orbit py on the 
collision manifold with < < \/2b there exist a manifold of orbits, lying 
on a cylinder, which eject from py. Similarly it can be shown that for every 
orbit Py, with —\/2b < v < 0, there exists a manifold of orbits, lying on a 
cylinder, which tend to py. 

If f = both types of manifolds exist, so po has a homoclinic manifold. 
Indeed, the equations that describe the manifold can be found explicitly: 
they have u = ±^/2b. With the energy relation we get 



V = ±V2r^h + 2r, (19) 
and using the equation of motion we obtain 



= ±rV2r^h + 2r. (20) 
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By integrating equation ( PP]) it is easy to find that 



4(r - To) 



(21) 



and 




'r 



2|/i| + (r-ro)2- 



2(r - To) 



(22) 



Furthermore 



f/(r - To) = ±V2b = u and ^{r - tq, ^o) = ©(^ - r^) - 9o, (23) 

where 0(r — tq) = a;(r — tq). As tq and 6*0 vary, equations (|2ll - p3D describe the 
entire 2-dimensional homochnic manifold. An orbit lying on the homoclinic 
manifold is represented in Fig. 3. Such an orbit is obtained by choosing 
^0 = 0; it ejects from the equator of the collision manifold, spiraling around it 
and moving upwards, then changes directions, goes downwards and upwards 
again, spiraling towards the periodic orbit Pq. 



Figure 3: An homoclinic orbit to po lying on the homoclinic manifold. This 
orbit spirals out of the equator of the collision manifold and then spirals back 
to it. 

The homoclinic manifold plays an important role in following section and is 
necessary for developing the generalization of the Melnikov technique. 
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5 A Generalized Melnikov Method 



Let X = {^{'T'):^{'^)^'^{'^)^U{t)) be the homoclinic orbit selected when we 
choose To = and 9o = 0. Consider solutions of the form 



r(r, To) = R{t - To) + r(r, tq) 
v{t, To) = V{t - To) + v{t, To) _ 

e{T, To, 60) = e(r - To) -60 + eir, ro) 

u{t,To) = U{t - To) + u{t,To). 
Let z = (f, V, 6, u), then the variational equation is 



z' = A{t)z + b(z, X, r, To, 9o, e) 



where 



A{t) 



and 



b(z, X, T, To, 9o, e) 



h 



( V i? \ 

l + ARh 

1 

/ 



r(r - To)v[T - To) 

-e((^cos2(e-^o + ^ 


y e(^ + 26) sin2(e-^o + 



The general solution of the variational equation 

i = $(t) ! ^~\s)hds. 



IS 



to 



(24) 



(25) 



(26) 



(27) 



see 



IB|), where $ is the fundamental matrix. If we let 



previous equation becomes 



Zi{t) = ^ij / Cj{s) ds 



(28) 
^-^h, the 

(29) 



to 



where cj = det Dj{t)/{det^){t) and Dj is the matrix obtained replacing the 
j-th column of $ with b. Furthermore the following formula for the trace 
holds: 

det$(r) =CeA^^^W'^^ (30) 
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One solution of the homogeneous part of the variational equation is given by 



X 

where 



[r - TO, Oo) = {R'{t - To), V'{t - To), 0'(r - To), U'{t - To)), (31) 



R' 
V 

&' 
U' 



4(r-ro) 



{2\h\ + {r 
2 



-ro)2)2 



4(T-ro)2 
(2|h| + (r-ro)2)2 



(32) 



2|/i|+(t-to)2 

±V2b. 
0. 

It is easy to check that other two independent solutions are (0,0,1,0) and 
(0,0,0,1). Knowing three independent solutions of a linear system, it is 
possible to find a fourth independent solution ip. This is achieved through 
the following lemma, which will be used to estimate how fast ip diverges. 

Lemma 1 Let z' = Az be the homogeneous part of ^B^). Given the three 
independent solutions above, a fourth is defined by 







ARh)/V'zA 
1 + ARh)zi 



R'{l+4:Rh)\~ 
— V' )^'^ 



^4 



R'zi + z/^ 
V'zi 

^3- 



(33) 



Proof: To find the fourth independent solution, we can use the "reduction to 
a smaller system," (see [|l^), whose direct application completes the proof. 
In particular it is useful to remark that we can always choose il)^ = 0, since 
{ipi,ip2, "^3, 0) is a solution that is independent from the other ones. To obtain 
necessary and sufficient conditions such that the negatively and positively 
asymptotic sets intersect transversely, we first obtain conditions for the exis- 
tence of solutions bounded on M for the non- homogeneous linear variational 
equation around x- 

For this, let B(R) = {h : R ^ R x R x x R bounded, continuous} with 
||b|| = sup,-g]g||b(r)|| for b G B(R). Then we have the following version of the 
Fredholm alternative for solutions bounded on M (see |^, Q for a similar 
approach) . 



Lemma 2 Let z G 



ixS^ xR and assume that z = m the expression 



of the function h. Then the variational equation 

A{t)z + b(z, X, r. To, 6*0, e) 



(34) 
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has a bounded solution if and only if 

/+00 
e- A ^'^^(^^'^^ R'{t - To) hix, r, To, ^o, e) dr = 0. (35) 
-oo 

The solution is unique and continuous and has the form z = L(b) + w, where 
L is a bounded linear operator, w = {0,0,9{to),u{oo)), when f(ro)-R'(ro) + 
v{to)V'{tq)) = 0, and satisfies the relation below, 

/+00 
h{x,'r^ro,eo,e) dr = 0. (36) 
-oo 

Proof: Using Lemma 1 it is easy to determine the behavior of as r ^ ±oo, 
precisely, 

5i ~ f "^Ai ~ T 

. ^2 ~ I ^2 ~ , . 

T ±00 < . , T ^ ±00 < , 4 (37) 

' ~ const. \ 1P3 ^ T 

T, V^4 ~ const. 

Using (^) and (0), the general solution of the complete (non-homogeneous) 
equation (P5|) can be written in integral form as 



f = R'i^A- J e- -^4 ^'■^('')'^''(7/;i62 - V^2&i) ds 

+ 4ji(^B + j\-^^^o^'^^^^''\B!h2 - V%) d^ 
v = V'(^A- j'^^- A^'-^('')'^''(V;i62 - ^2hi) ds^ 
+ V^2 (5 + j\- ^'J'^^''^''\B!h2 - V\) d^ 
e = ±V2b (^A - j\-^^o^'^'^^^''\^,b2 - ^2&i) ds 
+ ^3(^3 + ^%-^4^^^(^)'^''(i?'62 - V'b,) ds^ 

i-V'tPs ± V2bi^2)bi + (i?V3 ± V2biJi)b2 



(38) 



+ C- e~-''o 

J TO 

u = D + bi ds, 

J To 



TO 
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where, for notational convenience, we failed to mention the dependence on 
z, X, To, etc. 

Consider now the linearization of the problem ( pH]) around the solution 
z(r) = 0; in particular this amounts to deleting the high-order terms in 
the expression of b (i.e. bi = 0, etc.). Taking also into account the different 
behavior of the different solutions given in Lemma 1, it is easy to see that to 
have bounded solutions we need to require that 

^^(a- re^^4^^^('')">i62(x,s,ro,^o,e) rfs^ for z = 1, . . . , 4, (39) 



remains bounded as r ^ ±oo. More precisely z is bounded on [tq, oo) if and 
only if 

A= e-^=o^^^('')'^''^/,i62 ds (40) 

J TO 

and bounded on (— oo, Tq] if and only if 

/TO 
^-!:jrA{n)dn^^^^^^_ (41) 
■oo 

We also require 

u{±oo) = lim u{t) = D + hm / 64(x, s, Tq, 6'o, e) (is, (42) 

J TO 

where, obviously, n(oo) = n(— oo). The latter condition is not needed for the 
boundedness of the solution, but its role will be clear later when analyzing 
some properties of the negatively and positively asymptotic sets. It is easy to 
see that the above conditions are simultaneously satisfied both at r = — oo 
and at r = +oo if for some tq the following Melnikov-type conditions: 



+ 00 



/;^'-^(^)^>,62(x,s,ro.^o,e) ds = Q 



.Too (43) 
/ 64(x,s,ro,6'o,e) = 

J — oo 
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are fulfilled. Thus we can rewrite the general solution using (^) and, 
by neglecting to mention the dependence on x, s, etc., we obtain 



f = -R' Te-^^^^^^^^-^^ifea ds + ^1 (5 + Te-^^^^^^^^-^^i^'fes ds 

J oo \ J TO 

i = -V re-^4'^'^^('')'''7/;i62 ds + ^2 + f e' ^^^0^''^^^^"'' ds 

J oo \ J TQ , 

•J Tn 



u = u(oo) + / 64 (is 



(44) 



To obtain r(ro)-R'(ro) + ^(to) V(ro)) = we must have 



Moreover we also get 

C = 9{to) ±V2b e- -^4 ^''^^''^^''^ib2 ds (46) 
J 00 

and 

/•TO 

D = u{oo) + h ds. (47) 



00 



This uniquely defines B, C — 6'(ro), and D — u{oo) as continuous linear 
functionals on i3(]R). From ( PD we observe that the corresponding solution 
is of the form z = L(b) +w, where L is a bounded linear operator. It follows 
that this operator is continuous and hence the solution z = L(b) + w is 
continuous on i3(M). This completes the proof. 

To obtain necessary and sufficient conditions that the negatively and posi- 
tively asymptotic sets intersect, let us first consider all the solutions of (p5D 
which are bounded as r ^ — 00 and such that their angles remain close to 
the ones on the periodic orbit. The solution z is given by ( |38| ) satisfying (^TD 
and (|42|) with negative sign. In particular the solutions of the variational 
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equation that are bounded as r — — oo (i.e. which remain in a sufficiently 
small neighborhood of the periodic orbit as r ^ — oo) and with perturbed 
angles that do not drift but remain near the angles on the periodic orbit, 
must be on the negatively asymptotic set. In the same way, we obtain the 
positively invariant set from the solutions that remain bounded as r ^ oo 
and whose angles stay close to the one of the periodic orbit, which was in 
fact the reason why we required that condition (^) be satisfied. 
Moreover it is important to remark that the solution we found are not only 
bounded but also such that f— >-0, i)— >-Oasr— >oo and this is impor- 
tant since, on the collision manifold we have many periodic orbit and this 
condition is needed to show that the orbits are actually asymptotic to the 
equator. 

With the preparations above, we can now prove the following result. 



Theorem 3 System ^Tq ) has transversal homoclinic solutions if and only if 

dMi dM2 dMi dM2 



there exist Tq and a 9q such that 



TO=T( 




where 

/ + 00 
■oo 

/ + 00 
64(z*,r,ro,6'o,e) dr, 
-oo 



-S ^ 0, (48) 



(49) 



and z* is a solution ofz = L(b(z, r, tq, 6'o, e)) + w. Moreover if the perturba- 
tion is periodic we get infinitely many intersections. 

Proof: The stable and unstable manifolds intersect if and only if the solution 
([38| ) satisfies the Melnikov-like conditions (|35|) and (^2]) of Lemma 2. This 
was already proved in the case when b did not implicitly depend depend on 
z. But because of this implicit dependence we need to apply the implicit 
function theorem, which states that given z = L(b(z, r, Tq, 6*0, e)) + w with 
z — w = L(0, r. To, 6^0, 0) = 0, there exist a S and a unique solution z*(e, tq, 6^0) 
(that has continuous derivatives up to order 2 in tq, 6q, e) such that e < 0, 
|z| < 5 if the linearized operator z = L(b(0, r, tq, 6q, e)) + w is invertible. But 
Lemma 2 proved that such an operator is invertible. Moreover the homoclinic 
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solutions are transversal if and only if the integrals ( P^D have simple zeroes, 
as functions of Tq and 9q (see P, This concludes the proof. 

Unfortunately the Melnikov integrals of Theorem 3 are difficult to compute 
explicitly. To overcome this difficulty we need to rewrite these integrals to 
the first order approximation in e. Hence if we let z* = ex and b = ed with 
d = (c?!, ci2, (^3, (^4), the next result follows immediately. 



Corollary 1 System ( ffq j has transversal homoclinic solutions if and only if 

dMi dM2 dMi dM2 



there exist Tq and a 6q such that 



M,{T*,e*)=M2iT;,e;) = and 
where 



+00 

rTrA{s)ds 



r,=r* ^ 0, (50) 
80=80 



M,M) = / e-J^o^^-^^^"^ R'{r-ro) hixir - To),Qir - To) - 60) dr, 
J —00 
r+00 

M2{to, Oo) = / h{x{T - ro), 0(r - tq) - ^0) dr. 



(51) 

Moreover if the perturbation is periodic we get infinitely many intersections. 

Corollary 1 generalizes the Melnikov integrals obtained in [|1^, to non- 
hyperbolic whiskered tori (periodic orbits) in non-Hamiltonian systems. We 
remark that the second integral in (^1]) converges only conditionally. This is 
not a new feature of this non-Hamiltonian system since the same nuisance 
was present in []TU|, EBI . However some authors, more recently, found a way to 



write the Melnikov conditions for hyperbolic whiskered tori in Hamiltonian 
systems using only convergent integrals see |T^, It would be interesting 
to generalize those results to nonhyperbolic tori in non-Hamiltonian systems 
and to apply the newly developed technique to the problem under discussion 
in this paper. But this is not a project we aim to develop here. 
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6 The Melnikov Integrals 



Now we would like to apply Corollary 1 to our problem. The Melnikov 
conditions take the form 



Mi(ro,^o 



+00 



X cos^(ci;(r — tq) — 6*0) 



(52) 



dT = 



and 

M2(ro,^o) 



{R{t - To) + 2b) sin(2(cu(r - tq) - ^0)) dr = 0. (53) 



Let 9o = —6*0 — ujtq. With this assumption we can rewrite the first Melnikov 
condition as 

Ml = cos^ 9oI^ + sin^ Ool^ - sin 2^0/1", (54) 



where 



' I-= r+°°e-^^-;^W'^^i?/2'cos2ajr rfr 

J- J —00 



-00 
^00 



II = r +°° e--^^^o ^(^^'^^i^i?' sin^ UT dr 

J- J— 00 



-00 

hoo 



(55) 



jc _ r+co ^ 2 Jro ^^'^^'^'^RR' sin LOT COS LUT dr. 

V. J- J —00 

The second Melnikov condition can be expressed as 



where 



M2 = cos 2^0/2 + sin 260!^, 



= i J^^{R + 2b) sin 2ujT dr 
II = 1 J^^{R + 2b) cos 2uT dr. 



(56) 



(57) 



All the integrals above can be computed using the method of residues. 
Straightforward computations give 



and 



-F - — 

' " \h\ 



(r-ro)cos2cjr , vr sin(2curo)e- W^l'*! 
(2|/z| + (r-ro)2)2 



(r — To) sin CJT cos cur 
(2|/.| + (r-ro)2)2 



dr 



\h\^/2\h\ 
7rcos(2a;ro)e-^"^ 



(5J 



(59) 
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Particular care is needed when integrating and since they converge only 
conditionally. To obtain computational convergence, we choose the limits in 
1 2 such that 



J^' = lim / ib + — — : -77 1 sin 2iUT dr 



2 N^^J_^^/2cu\ 2|/i| + (r-ro)2 
7rsin(2cjro)e-2^v^ 



(60) 



The integral was also computed using the method of residues. Similarly, for 



I2, we have 



Ntt/2u 



/r, = lim / \h + -rr^, -, I cos 2uJT dr 

~Nn/2uj 



7rcos(2curo)e-2-V^ 



(61) 



and thus 



Ml = M2 = sin(2(cjro + ^0)) rT^- (62) 

a/2 /i 



We therefore have only one independent condition; this is clearly a conse- 
quence of the energy relation. 

We can find simple zeroes when sin(2(a;ro + 6*0)) = 0, i.e., for —{lutq + ^o) = 
^0 = ±W2 for A; = 0,1,2, .... 

Hence, by Corollary 1, we have proved the existence of an infinite sequence of 
intersections on the Poincare section of the negatively and positively asymp- 
totic sets of the periodic orbit and the existence of homoclinic orbits leaving 
the equator of the collision manifold and going back to it. This situation 
is clearly reminiscent of the chaotic dynamics described by the Poincare- 
Birkhoff-Smale theorem in terms of symbolic dynamics and the Smale horse- 
shoe. Unfortunately this theorem cannot be directly applied, nor can the the- 
orems proved in |]l|, since the Poincare-Birkhoff-Smale theorem considers hy- 
perbolic fixed points while the arguments in apply to area-preserving dif- 
feomorphisms. However the arguments contained in those theorems strongly 
suggest the occurrence of a chaotic dynamics. 

Moreover it is easy to verify, and interesting to remark, that the orbits we 
found above are not S'o-symmetric, where the So symmetry is defined by 
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So{r,v,9,u,T) = {r, —V, —u, —t) (see |]TU[) and an orbit 7(r) is said to be 
S'o-symmetric if 5*0 (7 (r)) = 7(r). Indeed an orbit is So-symmetric if and 
only if it has a point on the zero velocity curve, i.e., if there is a r such that 
v(t) = u{t) = (see ||2^). But this cannot happen in our problem because 
the unperturbed solution verifies u = ±\/2b. Thus for e small enough the 
perturbed orbit can never have u = 0. 

We can now summarize the above discussion as follows: 

Theorem 4 Let us consider the anisotropic Manev problem given by the 
equation of motion ^ with the energy relation jQ). Then there is an infi- 
nite sequence of intersections in the Poincare section of the negatively and 
positively asymptotic sets of the periodic orbits at the equator of the collision 
manifold (possibly giving rise to a chaotic dynamics) . Furthermore there ex- 
ist the homoclinic non So-symmetric orbits to the periodic orbit described 
above. 



7 Periodic Solutions 

We now return to the original Cartesian coordinates, which are more conve- 
nient for the purpose of finding certain periodic solutions. Let us first notice 
that the equations (|^) admit the following symmetries: 

So{x,y,Pa:,Py,t) = {x,y,-p^,-py,-t), 
Si{x,y,p^,Py,t) = {x,-y,-p^,py,-t), 
S2{x,y,p^,Py,t) = {-x,y,p^,-py,-t), 
S3{x,y,p^,Py,t) = {-x,-y,-p^,-py,t), 
Si{x,y,p^,Py,t) = {-x,y,-p^,py,t), 
S5{x,y,p^,Py,t) = {x,-y,p^,-py,t), 
S&{,x,y,p^,Py,t) = {-x,-y,p^,py,-t), 

which are the elements of an Abelian group of order eight, isomorphic to Z2 x 
Z2 X Z2, that is generated by 5*0, 5*1, S2 (see (The symmetry 5'o is the one 

denoted by 5*0 in the McGehee coordinates of the previous section.) To obtain 
certain families of periodic solutions, we will use the symmetries 50,51 and 
5*2 in connection with the variational principle according to which extremum 
values of the action integral yield periodic solutions of the equations (0). To 
reach this goal we first need to introduce some notations. 
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Let C°°([0,T],M2) be the space of T-periodic cycles / : [0,T] ^ R'^. 
Define the inner products 

(63) 

{f,9)m = {f,9)L^ + {f,9)L^, 

and let || -111,2, || ■ be the corresponding norms. Then the completion 
of C°°{[0,T],M?) with respect to the norm || ■ 11^2 is denoted by and it 
is the space of square integrable functions. The completion with respect to 
II ■ II is denoted by and is the Sobolev space of all absolutely continuous 
T-periodic paths that have derivatives defined almost everywhere (see 



Let Si([0, T], M^) denote the subset of formed by the S'j-symmetric paths, 
with i e {0, 1, 2, 3, 4, 5, 6}. It is easy to see that each is a subspace of H^; 
in fact they are Sobolev spaces and have many interesting properties. In the 
following we will restrict our attention to the spaces Sj with i = 0,1,2,6. 
Let us now prove the following result. 

Lemma 3 Let be defined as above, then the subspaces Sj of Si-symmetric 
paths with i = 0, 1, 2, 6 are closed, weakly closed, and complete with respect 
to the norm \\ ■ \\h^, and are therefore Sobolev spaces. Moreover 

i/^ = Si © S2 = So © Sg. (64) 

Proof: We first show an interesting fact: we can write / = {f^^f"^) as the 
sum of an 5*1 and an S'2-symmetric path. Indeed it is well known that we 
can write /i and /2 as the sum of an even and an odd absolutely continuous 
function, i.e. as fi = ft + f" and /2 = /I + /I- Using this idea we can 
write the path f(t) as the sum of an S'l-symmetric function, fs^ = (/i,/2)) 
and an 5'2-symmetric one, /sj = (/i, /!)• Now fix an element / G Si. Then 
{f^9)H^ = for every g G S2. This is because 

{f,9)m= {fi9i + f292) dt + / Ui9'i + 1292) dt, 
Jo Jo 

where the first integrand is an odd function and the second is an odd function 
almost everywhere. Thus the above scalar product is zero for every g G S2. 
Let us denote the space orthogonal to Si by Sf = {(7 G Si : {f,g)m = 
for every (7 G Si}. It is easy to see that S^*- is closed and that 5*2 C S^. 
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Now we need to show that 5*2 D S]*-. Assume there is /i G S]*- such that 
h and /i G S2. Then write h = hs^ +^S2 consider (/is^, hs^ +hs2)m, 
which means that (hs^ihs^) m = II^Sill = 5 > 0. But this contradicts the 
hypothesis that h G S^. Therefore E2 = S^. So S2 and consequently Ei are 
closed and such that = Si © S2. Moreover, since is a metric space, 
Si and S2 are complete. Also Si and S2 are weakly closed since they are 
norm-closed subspaces. The statements for Sq and Sg can be proved in a 
similar way. This completes the proof. 

Let us now introduce some new definitions. We will say that a path in Sj is 
of class Ln, n = 0, ±1, ±2, ±3, . . . , if its winding number about the origin of 
the coordinate system is n (i.e. if it makes n loops around the origin). The 
sign of n is positive for a counterclockwise rotation and negative otherwise. 
Consider the sets Si([0, T], M^yjo}). Notice that they are open submanifolds 
of the spaces Sj([0,T],M^) and that the family {Ln)nez provides a partition 
of those spaces into homotopy classes, also called components. Two periodic 
orbits of the isotropic Manev problem [n = 1), one of class Lg and the other 
of class L_9, are depicted in Fig. 4. 




(a) (h) 



Figure 4: Periodic orbits of the Manev problem: (a) 5'i-symmetric periodic 
orbit of class Lg, (b) S'2-symmetric periodic orbit of class L.g. Note that 
none of these two is S3 symmetric. 

The Lagrangian L(q, q) = T(q) + f/(q) of the anisotropic Manev problem 
given by system @ has the expression 

L{x,y,x,y) = l{x^ + f) + ^ 5 + ^TT — 2^ (6^) 
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and the action integral along a path / from time to time T, whose Euclidean 
coordinate representation is q = q(t) = {x(t),y{t)), takes the form 

AT{f )= r L{^{t)A{t))dt. 



Jo 

According to Hamilton's principle, the extremals of the functional At are 
solutions of the equations (0). Hence we want to obtain periodic solutions 
of @ by finding extremals of the functional A. For this we will use a direct 
method of the calculus of variation, namely the lower-semicontinuity method 
(see [^). In preparation of a satisfactory theory of existence, the notion of 
admissible function has to be relaxed since the extremals we obtain belong to 
a Sobolev space. Therefore the above method provides only "weak" solutions 
of our problem. To show that the paths are regular enough to be classical 



solutions, we need the following result, proved in 



Lemma 4 The critical points o/ At|sj([o t] r2\{o}) '^'^^ T -periodic solutions of 
equations 

In particular it is well known that if / is a minimizer of the action At in 
the space -ff^([0, t], M^) and if / has no collisions, then / is a T-periodic 
solution to (I). Collision have to be excluded because equations (H) break 
down at collisions and because the action is not differentiable at paths with 
collisions. In this paper we are interested to restrict ourself to the spaces Sj 
of S'j-symmetric paths for z = 0, 1, 2. The paths that are S'e-symmetric have 
to be excluded in the study of periodic orbits since Se-symmetric paths must 
intersect the origin and therefore encounter collisions. 

Now it is not obvious that a coUisionless minimizer in Sj is a periodic solution 
of system @j. However, according to the principle of "symmetric criticality" 
(see for example ^) this is actually true. Indeed, it can be proved that 
if / is a collision free path with dAt{f){h) = for every f & "Ei, then 
dATif)ih) = for all / e //^([O, T], M^) and thus / is a critical point in the 
bigger loop space (see 0). 

The only obstacle left for applying the direct method is the "noncompact- 
ness" of the configuration space. Indeed we want to exclude the possibility 
that the minimizer is obtained when the bodies are at infinite distance from 
each other or are collision paths. The first problem is solved restricting our- 
selves to non-simple cycles, i.e., to cycles that are not homotopic to a point 
and thus are not in the homotopy class Lq. The second problem is solved by 
the following result. 
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Lemma 5 Any family T of non-simple homotopic cycles in'Ei{[0,T],'R'^\{0}) 
fort = 0,1,2 on which J{f) = ^\q{t)\'^dt and E{f) = U{q{t))dt are 
bounded, is bounded away from the origin. 

The proof of this resuh follows from [|l^ if we remark that the anisotropic 
Manev potential is "strong" according to Gordon's definition and that the 
Lagrangian is positive. 

To apply the direct method we still need to recall some properties of lower 
semi continuous (l.s.c) functions. Let JF : X — > M be a real valued function on 
a topological space X. Then JF is l.s.c. if and only if JF~^(— oo,a] is closed 
for every a G M, in which case JF is bounded below and attains its infimum 
on every compact subset of X. Moreover when X is Hausdorff then compact 
sets are necessarily closed and thus we have the following result. 

Proposition 2 Suppose J-" : X W is a real valued function on an Haus- 
dorff space X and 

J-'~^{—oo,b] is compact for every real b. 

Then T is l.s.c, bounded below, and attains its infimum value on X. 



We can now prove the main result of this section. 



Theorem 5 For any T > and any n = ±1, ±2, ±3, . . . , there is at least 
one Si-symmetric (i = 0,1,2) periodic orbit of the anisotropic Manev problem 
that has period T and winding number n (i.e., belongs to the homotopy class 

Ln). 

Proof: Let X be a component of Sj([0, T], ]R^\{0}) for z = 0, 1, 2, that consist 
of non-simple cycles. Endow X with the weak topology it inherits from 
Ej ( [0, T] , M^) . Then X is a subset of an Hilbert space and it is weakly compact 
if and only if it is weakly closed. 

We wish to apply Proposition 2 with JF = At and thus we have to show that 
X n ATp^ {—oo, h] is a bounded and weak-closed subset of Sj([0, T], M^). 
Since J = Ax — E and > 0, we have E > and therefore 

J<h on AT~\-^,h]=AT~^[0,h], 

E = At - J <h on ^^^^[0,6]. ^ 
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Since J <b the elements of X are bounded in arc length, and from Lemma 5 
it follows that the elements of X are bounded away from the origin. Moreover 
the elements of X are non-simple and thus bounded in the C° norm and hence 
in the norm. This last fact combined with J < b shows that X is bounded 
in the || ■ norm. Thus also X fl At~^ {—oo, b] is bounded in the norm. 
Now suppose that {/„} = {(/^, /^)} is a sequence in X DAt^^IO, b] that con- 
verges weakly to a cycle / G Zlj([0, T], M^) for i = 0, 1, 2. From general princi- 
ples, ll/nll/fi is bounded and ||/rt||L2 — > ll/nlU^ because weak Ej- convergence 
implies C°-convergence. Since J{fn) = l/2||/n||//i — l/2||/n||^2 it means that 
J{fn) is bounded and since E < b on At~^[0, 6] it follows that {£"(/„)} is 
bounded. Moreover, Lemma 5 guarantees that the functions /„ are bounded 
away from the origin so that / is homotopic to the /„ in ]R^\{0}. Therefore 

fex. 

To complete the proof we have to show that / G At ^[0,6]. We know that 
E{fn) E{f) since weak convergence in Sj implies C°-convergence. For 
each n let 

1 1 

9n{t) = . ........ + 



and denote 

1 1 



9{t) 



Each Qn is of class since Aj-{fn) < oo. This implies that the set of all t 
for which fn{t) = has zero measure, otherwise the integral of gn{t) would 
be unbounded. So (7„(t) — ^ g{t) almost everywhere. Also Quit) dt < 
Arifn) ^ b. By Fatou's lemma it follows that g is and that 

g(t) dt = liminf 5f„(t) dt < liminf / g-nit) dt. 
Jo Jo 

Now we can use the fact that the norm is weakly sequentially lower semicon- 
tinuous (see p5i), thus 



= WfWm - WfWh < liminf WfnWl. - ml. = liminf 
where the last equality holds since {/„} converges strongly to / in L^. Con- 
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sequently 

^T(/)=^ii/iii.+ r9{t)dt 

1 fT ^6^) 

< liminf -||/n||i2 + liminf / g^it) dt < liminf y4r(/n) < b. 
^ Jo 

Relation (|67|) now implies that / G y4T'^[0,6]. This completes the proof. 

Recall now that two intersections of every Si-symmetric (S'2-symmetric) or- 
bit with the X axis [y axis) must be orthogonal. To distinguish them from 
accidental orthogonal intersections, which do not follow because of the sym- 
metry, we will call them essential orthogonal intersections. From the proof 
of Theorem 5 and obvious index theory considerations, the following result 
follows (see also Fig. 4). 



Corollary 2 If the essential orthogonal intersections with the x-axis (y-axis) 
of an Si-symmetric (S2-symmetric) periodic orbit lie on the same side of the 
axis with respect to the origin of the coordinate system, then the orbit has 
an even winding number. If the essential orthogonal intersections are on 
opposite sides with respect to the origin, then the periodic orbit has an odd 
winding number. 

Since the symmetries 5*0, 5*1 and S2 generate the entire symmetry group, it 
is clear that Theorem 5 captures all periodic orbits with symmetries. This 
result, however, does not tell if other symmetric periodic orbits exist beyond 
the the ones with 5*0, Si and 5*2 symmetries. Let us therefore end our paper by 
proving that S'3-periodic orbits do indeed exist. In fact they form a rich set if 
compared to the one of S'3-symmetric orbits of the anisotropic Kepler problem 
(given by (|l]) with 6 = 0), which contains only circular orbits. We will show 
that in our case each homotopy class n = Ak + 1, k integer, contains 
at least one S'3-symmetric periodic orbit. Other homotopy classes may have 
S's-symmetric periodic orbits, but our approach proves their existence only 
for winding numbers of the form n = 4k + 1, k integer. 
We consider the set of all paths with one end on the x axis and the other on 
the y axis of the coordinate system. As in the case of periodic cycles discussed 
in the first part of this section, for a given T' = T/4 > this set can be 
endowed with a Hilbert space structure, the completion of which is a Sobolev 
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space. We further divide this space in homotopy classes n = 0, ±1, ±2, . . . 
according to the winding number n. 

Using the boundary conditions, it is easy to see that in each class the 
minimizer of the action is a an arc orthogonal to the x and y axes. Its ex- 
istence and the fact that it is a solution in the classical sense can be proved 
in a similar way as we did for periodic cycles. Once obtaining such a solu- 
tion with ends on the x and y axes, we can use the 5'3-symmetry and the 
orthogonality with the axes to complete this solution arc to a periodic orbit 
of period T > 0. The symmetry implies that the winding number is of the 
form n = Ak + 1, k integer. This is because if, for example, a solution arc 
with the ends on the x and y axes has a loop around the origin, then the 
corresponding periodic orbit has four loops around the origin. We have thus 
obtained the following result. 

Theorem 6 For any T > and any n = Ak + 1, k integer, there is at least 
one S^-symmetric periodic orbit of the anisotropic Manev problem that has 
period T and winding number n (i.e., belongs to the homotopy class L^). 

It is interesting to note in conclusion that if viewing the anisotropy parameter 
as a perturbation and the anisotropic Manev problem as a perturbation of the 
isotropic case (see Section 4), the S'j-symmetric {i = 0,1,2) periodic orbits 
of the isotropic problem are deformed but not destroyed by introducing the 
anisotropy, no matter how large its size. This shows that the Si symmetries 
{i = 0, 1, 2) play an important role in understanding the system and are an 
indicator of its robustness relative to perturbations. 
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